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Abstract 

We re-examine the big bang nucleosynthesis (BBN) bounds on the mixing 
of neutrinos with sterile species. These bounds depend on the assumption 
that the relic neutrino asymmetry L v is very small. We show that for L v 
large enough (greater than about 10~ 5 ) the standard BBN bounds do not 
apply. We apply this result to the sterile neutrino solution to the atmospheric 
neutrino anomaly and show that for L v > 7 x 10~ 5 it is consistent with BBN. 
The BBN bounds on sterile neutrinos mixing with electron neutrinos can also 
be weakened considerably. 



The solar neutrino deficit Jl[, atmospheric neutrino anomaly 0, and 
LSND experiment |3|] can all be viewed as evidence for non-zero neutrino 
masses and oscillations. It does not seem possible to explain all these anoma- 
lies with the three known neutrino species and thus new neutrino species 
might exist. Given that new ordinary weakly interacting neutrino species 
are ruled out by LEP, sterile neutrinos (i/ s ) are a natural candidate. There 
are essentially two types of sterile neutrinos that can be envisaged. Firstly, 
there are sterile states which either have no gauge interactions, or interactions 
which are much weaker than the usual weak interactions [|J. Alternatively, 
it is possible to envisage neutrinos which do not have significant interactions 
with ordinary matter but do have significant interactions with themselves. 
An interesting example of the latter is given by mirror neutrinos which inter- 
act with themselves only through mirror weak interactions which have the 
same strength as ordinary weak interactions PJ. 

However, for both sterile and mirror neutrinos there are apparently quite 
stringent bounds if they are required to be consistent with standard big bang 
cosmology. Assuming that the number of effective neutrino species present 
during nucleosynthesis is bounded to be less than 4, then the mixing angle 
(#o) and the squared mass difference (5m 2 ) for a sterile neutrino mixing with 
one of the known neutrinos is bounded by (assuming 5m 2 > 0)0: 

5m 2 sin 4 29 ~ 5 x 1(T 6 eV 2 , v = u e , 

5m 2 sin 4 29 ~ 3 x KT 6 eV 2 , v = v^ T . (1) 

These bounds arise by demanding that oscillations do not bring the sterile 
neutrino into equilibrium with the known neutrinos. Electron neutrinos must 
also not be depleted too much by oscillations after decoupling, during the 
BBN epoch, because then the freeze out temperature for neutron-proton 
transitions is increased. For maximal mixing, the bound on 5m 2 is extended 
to 5m 2 ~ 1CT 8 eV 2 ||. These "bounds" would appear to exclude the region 
of parameter space required to explain the atmospheric neutrino anomaly 
in terms of — u s oscillation (5m 2 ~ 10~ 2 eV 2 , sin 2 29q ~ 1), and would 
restrict the parameter space required to explain the solar neutrino deficit in 
terms of v e — v s oscillation. An important assumption in deriving the bounds 
Eq.(l) is that the relic neutrino asymmetries could be neglected. However, 
the neutrino asymmetries cannot be measured, and at present the origin of 
particle asymmetries is not fully understood. Only the extremely weak bound 
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L Ua ~ 10 3 can be derived by demanding that the neutrinos do not violate 
the upper limit on the total energy density of the universe. The purpose of 
this letter is to re-examine the BBN bounds on ordinary-sterile neutrino mass 
and mixing for arbitary neutrino asymmetries. In particular, we will show 
that for neutrino asymmetries larger than about 7 x 10~ 5 the standard big 
bang model is consistent with sterile neutrinos mixing with muon neutrinos 
with parameters suggested by the atmospheric neutrino anomaly J7|. 

Let us first examine the ordinary neutrino (v a , a = e,/i,r) oscillating 
with a sterile neutrino (u s ) in vacuum. Oscillations can occur if the weak 
eigenstate neutrino and sterile neutrino are each linear combinations u a = 
cos OqUi + sin 9qv 2 and v s = — sin #0^1 + cos #0^2 of mass eigenstates z/ 12 . An 
ordinary neutrino of momentum p will then oscillate in vacuum after a time 
t with probability 

|(z/ a (t)|z/ s )| 2 = sin 2 2^ sin 2 f-^V (2) 



Lose 



where 



L - 2P = 1 (3) 

Losc ~ Sri 2 " ~ v { ) 

However, in the early universe oscillations occur in a plasma. For v a — v s (a = 
e, fi, t) oscillations in a plasma of temperature T, the matter and vacuum 
oscillation parameters are related by 

• 2 na sin 2 2^o 

sin 2u„ 



1 - 2z cos 29 + z 2 

A 2 m = Al(l-2zcos26 + z 2 ), (4) 

where z = 2{p)(V a — V s )/5m 2 and (V a ^ s ) are the effective potentials due to the 
interactions of the neutrinos with matter ({p) ~ 3.15T). For a truly sterile 
neutrino V s = 0. (For neutrinos which have only self interactions, e.g. mirror 
neutrinos 0, V s can be non-zero. We will comment on this case later.) For 
a weak eigenstate neutrino (z/ Q , a — e, fi, r), V a is given by (l^, [TTJ, [L2| 



V a = V2G F N, (lM - ^ j , (5) 
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where Gp is the Fermi coupling constant, My/ is the W boson mass, A a is 
a numerical factor given by A e ~ 55 and A^ T ~ 15.3 [10, 11|. The neutrino 
asymmetry is given by 

L (Q) =L a + L Ue +L^+L UrJ (6) 

where L a = (iV Q — Ng^/N^. Since we will be interested in the case where 
the asymmetries are large (of order 10 -5 or more) we have neglected the 
asymmetries in the electrons and protons/neutrons since these are known to 
be small (~ 1(T 10 ). 

For L (q) non negligible, the bounds of Eq.(l) can be weakened consider- 
ably. This is because the oscillation probability depends on through the 
dependence of sin 2 29 m on in Eqs.(4,5). The condition that the sterile 
neutrinos not come into equilibrium is that the interaction rate for sterile 
neutrinos is less than the expansion rate, i.e. T Us < H. We will assume 
that there are essentially no sterile neutrinos initially. The rate of produc- 
tion of sterile neutrinos is given by the interaction rate of ordinary neutrinos 
multiplied by the probability that the neutrino collapses to the sterile eigen- 
function, i.e. 

r u s {t) = (Pu a -> ^s)coll^u a , (7) 

where Y Va = y a G 2 F T 5 , (y e ~ 4.0 and y^ T ~ 2.9) f§ and (P Ua -> u s ) coll is 
given by 

(P Ua -> v s ) coU = sin 2 26> m (sin 2 — ^y), (8) 

Lose 

where x is the distance between collisions. Note that (x) = L int = l/Y Va 
where Li nt is the mean distance between interactions. From Eq.(4) it is easy 
to see that the production rate of sterile neutrinos is significantly suppressed 
(because sin 2 29 m sin 2 20 o ) for temperatures above about 12 MeV (given 
8m 2 ~ 10~ 2 eV 2 ), independently of the magnitude of (except in a reso- 
nance region where sin 2 29 m = 1, as we shall discuss later). Below about 12 
MeV, sin 2 26' m approaches its vacuum value (unless is non negligible). 
In the standard scenario^, it is in this region where oscillations can occur 
and potentially bring the sterile neutrino into equilibrium. 

However, for non negligible, the production rate of sterile neutrinos 
can continue to be suppressed in the region below 12 MeV. Assuming \L^ | > 
10~ 5 (which we will discuss later), the condition Li nt /L^} ^> 1 always holds 
off resonance (for temperatures above the decoupling temperature). This 
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means that (sin 2 (L inf /L^)) averages to 1/2. Using this result, and Eqs.(7,8) 
[with sin 2 29 m given by Eq.(4)], we can now calculate the production rate of 
sterile neutrinos for the general case with non-zero. Demanding that this 
interaction rate be less than the expansion rate, i.e. T Us ~ H ~ 5.5T 2 /Mp, 
we find for large (where the second term in Eq.(5) can be neglected) and 
maximal mixing 



2 / y a M P GpT 3 



(6m 2 ) 2 Z 79G 2 F T 2 N 2 [L™\ 1 ^ ^ - 1 1 , (9) 

where M P ~ 1.2 x 10 22 MeV is the Planck mass. Note that in the case 
where the mixing is not maximal, the bound is even more stringent. Clearly, 
oscillations from ordinary to sterile neutrinos which occur after the kinetic 
decoupling temperature (T dec ) do not significantly affect the energy density 
of the universe. For temperatures above T dec , the most stringent bound 
occurs at the decoupling temperature. For z/ e , T dec ~ 2.6 MeV and for z^ r , 
Tder. — 4.4 MeV, which leads to the bounds: 



v 



\5m 2 \ < 4 x 10 2 \L {e) \eV 2 , 

\5m 2 \ < 1.6 x 10 3 |i>' T) |e1/ 2 , v = v^ T . (10) 

These bounds replace those of Eq.(l) in the case where \L^ a ' | is large [\L^ | > 
lO" 5 ]. 

In the case of electron neutrino oscillations into sterile neutrinos a more 
stringent bound comes by requiring that the electron neutrinos are not de- 
pleted significantly down to temperatures where the protons and neutrons 
go out of equilibrium (which is about 0.8 MeV in the standard scenario). 
For maximal mixing, the bound is 5m 2 < 10~ 8 eV 2 if is negligible|J. 
However in the case where \L^\ is not negligible (> 10~ 5 ), the situation is 
changed. If we demand that sin 2 26 m ~ 1/10 for T > 0.8 MeV (which means 
that the v e — v s oscillations are severely suppressed for T > 0.8 MeV) we find 
in the case of maximal mixing that 

\5m 2 \ ~2 v / 2|L (e) |G F V 7 T ~ 3.3|L (e) | eV 2 . (11) 

This bound for v e — v s oscillations is clearly more stringent than Eq.(10). 

In the above analysis we have assumed that the value of L^ a ' is fixed. 
However in reality l/ a ) is in general not constant. Oscillations can change 
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its value. There are only two regions where oscillations are important. 
First, there is the region near the high temperature resonance (in this re- 
gion sin 2 29 m = 1). The other region where oscillations can significantly 
change is at the low temperature region T rs Tdec- Oscillations can 
change in this region because the production rate of sterile neutrinos is 
not so strongly suppressed (recall sin 2 29 m — ► sin 2 29 as T — > 0). We now 
consider each of these regions in turn. 

We will assume for definiteness that ife*' is positive (unless explicitly 
stated otherwise) . The production rate of sterile neutrinos is given by Eqs. (7,8) 
with sin 2 2# m given by Eqs. (4,5). The condition for v a — v s oscillation reso- 
nance {9 m = 7r/4) is that 

(V a - V s ) = A cos 26»o, (12) 



implying that the resonance temperature is[ 12 



T2 _ L^M 2 W A cos29 M 2 v 

A a A a V2G F N y ' 1 ' 

Observe that strictly, the right-hand side of Eq.(13) is a function of tempera- 
ture so that we must solve Eq.(13) for the resonance temperature. However, 
it turns out we will only be interested in the high temperature resonance and 
quite large values of L^ a \ i.e. ~ 10~ 5 , and in this case the second term 
on the right-hand side of Eq.(13) can be neglected (for 5m 2 < 1 eV 2 ). 

Observe that for > 0, 5m 2 > 0, there is no resonance for antineu- 
trinos while for L^ a ' > 0, Sm 2 < the resonance for antineutrinos occurs 
at very low temperatures. [Note that for 8m 2 < 10~ 2 eV 2 , L^ a ' > 10~ 5 this 
low temperature resonance occurs at temperatures below the kinetic decou- 
pling temperature and thus it can be neglected in our analysis.] Of more 
importance is the "high temperature" resonance which, for > only oc- 
curs for neutrinos (for L^ a > < the high temperature resonance only occurs 
for antineutrinos). The effect of the high temperature resonance is rather 
interesting. For an initial less than a certain "critical" value (to be 
determined later) U- a ' will evolve to zero. This is essentially because the 
oscillations near the resonance are so numerous as to continually lower the 
value of \L^\ and hence also the resonance temperature, so that the system 
cannot actually pass through the resonance. 
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However, for large enough there will be a critical point where the 
expansion of the universe is more important than the change in due 
to oscillations. This behaviour has been studied numerically in Ref . [O, O . 
Below we show how this behaviour can be understood and we derive an 
analytic approximation for the critical value of L^ a \ 

If the change in the resonance temperature due to oscillations is greater 
than the width of the resonance, then the resonance will be moved to lower 
and lower temperatures, until the lepton number is reduced to near zero. 
However if the change in the resonance temperature due to oscillations during 
the resonance is less than the width of the resonance [ j"4"| , then this will ensure 



the system passes briefly through the resonance. Having passed through the 
resonance, the value of will not change significantly until much lower 
temperatures, as discussed earlier. The condition that the system passes 
through the resonance is that 

5T res £ AT, (14) 

where 5T res is the change in the resonance temperature due to the oscillations 
as the system passes through the resonance, while AT is the width of the 
resonance. The resonance temperature T res is related to the lepton number 
asymmetry through Eq.(13) so that 

M w 5LW 

dl res = (15) 



Now, 5L^ is proportional to the reaction rate for ordinary neutrinos 
converting into sterile neutrinos multiplied by the time it takes for the system 
to pass through the resonance At, i.e. SL^ = — 3T Us At/4:. Note that the 
resonance width measures defined in terms of the temperature and the time, 
AT and At, are related to each other using the time-temperature relation of 
the early universe: t ~ Mp/llT 2 implies At ~ — MpAT/5.5T 3 . Hence, the 
condition that the system will pass through the resonance is that 

3M w M P T Us < i 



44T 3 y A a L^ 

Note that for > 1(T 5 , T res > 35 MeV. It is easy to verify that at the 
resonance, L int /L^ < 1 and hence sm 2 (L int / L^})) ~ L 2 nt /L^ 2 . Using 
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this result, and Eqs. (7,8), we can calculate T Ug at the resonance: 

A 2 

^ |res= y a G 2 F T^ (17) 
Substituting T Va \ res into Eq.(16), we obtain 
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^ -\Zxl&ylG%M%) ■ [l °> 

Note that L^ t is independent of the vacuum mixing angle 6*0. For < L^} t 
the change in the resonance temperature is greater than the width of the 
resonance. This means that the resonance dynamically evolves to later and 
later times, with moving closer and closer to zero. For > L^j t the 
change in the resonance temperature is less than the width of the resonance, 
so the system passes through the resonance. Having passed through the 
resonance the value of remains approximately unchanged (until much 
later times as will be discussed later). Thus requiring > L^/ it we find 

L (e) > 9 x KT 6 for 5m 2 < 10~ 4 eV 2 , 

i> r) > 1.8 x 1(T 5 for 5m 2 < KT 2 eV 2 . (19) 
It is useful to compare our analytic expression Eq.(18) with the numer- 



ical work of Ref.[13|. They find numerically that for 5m = 10 eV an 
initial asymmetry of L^ e ' = 10 -5 remains unchanged on passing through 
the resonance. They also obtain that for 5m 2 = 10~ 3 eV 2 an initial asym- 
metry L$ = 10 -5 leads to evolve to zero. These results are consis- 
tent with our analytic expression Eq.(18). [Putting in 5m 2 = 10 -4 eV 2 in 



Eq.(18), we find L^ t - 9x 1(T 6 < 1$ while for 5m 2 = 10~ 3 eV 2 we find 
4 e r ! t ^2.1xl(T 5 >4 e) .] 

As a consistency test, we should check that the change in L^ a ' on passing 
through the resonance is small compared with the initial value of L^ a \ To 
work out the change in L^ a ' on passing through the resonance, we must work 
out the width of the resonance. Note that the width of the resonance is larger 
than might be expected. Actually, the interaction rate at the centre of the 
resonance 8 m = it/4 is equal to the interaction rate anywhere in the region 
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where L)™} = A m x ^> L int . To see this observe that the interaction rate is 
given by, 

r - = r - sin ' 2 M » sm2 (£) = ^ T5 f 8m2 (dN ■ (20) 

Thus, r^ s = T Va \ res [defined in Eq.(17)] provided A m < = y a G 2 F T 5 . 
Using Eqs. (13-15), it is easy to show that the change in on passing 
through the resonance is related to the resonance width by = 
2AT/T. Calculating AT using A m (T res + AT/2) ~ y a G 2 F T r 5 esJ we find that 
SL^/L^ < 0.04 (0.10) for v a = u e (v a = v^ T ). Thus, the the change in 
through the resonance is always at least an order of magnitude smaller 
than Z>). 

Note that the calculation of L^ it assumes that the sterile neutrino has no 
significant interactions. In the case of a neutrino which has only significant 
interactions with itself (such as a mirror neutrino H) then V s is unequal to 
zero. In this case L^ a ' does not evolve to zero but evolves so that V a — V s ~ 0. 
This does not affect the resulting analysis, since the interaction rate depends 
on V a — V s rather than V a . 

We now consider the low temperature region T « T^ ec in which oscilla- 
tions can also potentially erase L^ a \ To calculate SL^ in this region, we 
must integrate <5L (q) = -3(r„. - T Va )6t/4, using T = (l/2)r„ a sin 2 O /(1 - 
2z cos 29 + z 2 ) where z differs between v and V by L^ a ' — > —L^ a \ Approxi- 
mating 1 — 2z cos 29 + z 2 by z 2 , which is approximately valid for Lr- a > > 10 -5 
and 5m 2 < 10~ 2 eV 2 (10~ 5 eV 2 ) for a = fi, r and T > 4.4 MeV (a = e and 
T > 0.8 MeV), we find that 

fr(«n 4 ^M«) \ 4 3 sin 2 e A a y a M P (5m 2 ) 2 5 4 

V^finalJ — ^initial J 22M 2 T 3 

where the last inequality comes by demanding that Lf*^ > 10~ 5 so that the 
bounds Eqs. (10, 11) remain valid. Note that the strongest tendency towards 
erasure of occurs at the low temperature end of the integration region. 
For a = (jl,t, Tf = T dec ~ 4.4 MeV. For a = e, we require Tf ~ 0.8 MeV 
because we require that not be erased above the temperature in which 
the protons and neutrons are kept in equilibrium. Evaluating Eq.(21) we 
find that in the most stringent case of maximal mixing: 

^initial > 4 x 10_5 for &m 2 < 10 " 4 eV ^ 
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L u£L > 7 x 10_5 for 5m2 < 10 " 2 ey2 - ( 22 ) 

These bounds are slightly more stringent (for maximal mixing) than those 
obtained earlier from requiring that not be erased due to the high tem- 
perature resonance. 

Thus, we conclude that if L^ a ' satisfies Eqs.(19) and (22) then L^ a > is not 
erased either at high temperatures or at low temperatures, and hence for 
satisfying Eqs.(19) and (22) the bounds Eqs.(10, 11) hold. We conclude 
that for ZA a ) satisfying Eqs.(19) and (22) the sterile neutrino solution to 
the atmospheric neutrino anomaly is consistent with BBN. Also, the large 
angle MSW v e — v s oscillation solution to the solar neutrino problem is also 
consistent with BBN as is the maximal mixing vacuum oscillation solution 
(see e.g. Ref.§). 
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